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FOREWORD 


This  sponsored  research  program  closely  accompanied  research  activities  of  the  DURINT 
grant  sponsored  by  the  U.S.  Army  Research  Office.  The  specific  research  objectives  of  this 
grant  were  on  mathematical  modeling  of  ultrafast  condensed  phase  reactions  with  and  without 
pressure  gradients  generated  by  the  reaction  and  on  investigation  of  combustion  characteristics 
of  reacting  systems  consisting  of  nanopowders.  One-dimensional  mathematical  models  were 
derived  and  solved  numerically.  These  models  predict  qualitatively  combustion  front 
propagation  in  such  reacting  systems.  It  was  also  found  that  surface  functionalization  of 
aluminum  nanopowders  significantly  slows  down  its  reaction  with  moisture.  It  was 
demonstrated  experimentally  that  coated  aluminum  nanopowders  with  oleic  acid  or  silanes  react 
even  better  with  copper,  molybdenum,  and  tungsten  oxides.  This  interesting  discovery  can  be 
explained  by  more  ultimate  mixing  of  heterogeneous  reactants.  In  addition,  results  on 
characterization  of  nanopowders  and  mixing  of  nanoreactants  are  presented  as  well. 


LIST  OF  FIGURES 

Figure  1.  TEM  image  of  passivated  aluminum  nanopowder. 

Figure  2.  X-ray  patterns  of passivated  aluminum  and  products  of  combustion  in  oxygen 
(alumina)  and  in  air  (alon). 

Figure  3.  Schematic  of  the  glassware  for  the  determination  of  the  reactive  aluminum 
content  by  the  volumetric  method. 

Figure  4.  Schematic  of  the  stopper  used  in  neck  6.  Sample  cup  is  a  metal  piece  bonded 
to  the  end  of  the  glass  stopper. 

Figure  5.  Schematic  of  an  aging  chamber. 

Figure  6.  Aging  results  for  uncoated  Al  (UFAl  1-65-38). 

Figure  7.  Aging  results  for  coated  and  uncoated  Al  at  43  %RH. 

Figure  8.  Aging  results  for  coated  and  uncoated  Al  at  75  %RH. 

Figure  9.  Aging  results  for  coated  and  uncoated  Al  at  84  %RH. 

Figure  10.  Aging  results  for  coated  and  uncoated  Al  at  97  %RH. 

Figure  11:  Results  for  aging  of  Nanotechnologies  and  Technanogy  aluminum  samples  in 
97  %RH  at  40  °C.  Lines  shown  are  trendlines  only  and  each  data  point  is  an 
average  of  two  volumetric  method  measuremen  ts. 
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Figure  12.  Al-TiO 2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl 
sulfosuccinate  as  surfactant  ( SE-image  10000  X). 

Figure  13.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl 
sulfosuccinate  as  surfactant  (BSE-image  10000  X). 

Figure  14.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (SE-image  10000  X). 

Figure  15.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (BSE-image  10000  X). 

Figure  16:  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl 
sulfosuccinate  as  surfactant  (Cameo-image  10000  X). 

Figure  17.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (Cameo-image  10000  X). 

Figure  18.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl 

sulfosuccinate  as  surfactant  (Elemental  Mapping  10000  X SE-BSE-Al  and  Ti 
mapping  image). 

Figure  19.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (Elemental  Mapping  10000  X SE-BSE-Al  and  Ti  mapping  image). 

Figure  20.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl 

sulfosuccinate  as  surfactant  (Three  line  scans  of  10  pm  at  10000  X). 

Figure  21.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (Three  line  scans  of  10  pm  at  10000  X). 

Figure  22.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl 

sulfosuccinate  as  surfactant.  Sample  after  three  line  scans  of  10  pm  at 
10000X. 

Figure  23.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl 

sulfosuccinate  as  surfactant.  Sample  after  three  line  scans  of  10  pm  at 
10000X.  The  X-  ray  spot  analyses  at  four  points  for  aluminum  and  titanium 
show  the  following  results  #1  ( Al :  44.7  Ti:  55.3),  #2  ( Al :  48.4  Ti:  51.6),  #3  (Al: 
48.8  Ti:  51.2),  #4  (Al:  40.8  Ti:59.2). 

Figure  24.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant.  Sample  after  the  three  line-scans  of  10  pm  at  10000  X. 


Figure  25.  Images  of  burning  loose  nanopowders  of Al-MoO  3  system  recorded,  with 
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perforated  baffles,  using  high-speed  camera. 

Figure  26.  Results  for  the  measurement  of propagation  velocity  of  Al/CuO  as  a  function  of 
weight  percent  coating  on  the  aluminum  nano-powder.  The  coatings  used  were 
Z6124  silane  and  oleic  acid.  Each  data  point  is  the  average  of  4  measurements. 

Figure  27.  Propagation  velocity  as  a  function  of  coating  concentration.  Coatings  were  added  at 
the  time  of  mixing  of  Al/CuO  reactants 

Figure  28.  Al/CuO  ignition  delay  results  for  coated  aluminum  sample.  Coatings  studied  were 
Z6124  Silane  and  oleic  acid. 

Figure  29.  Schematic  of  the  cylindrical  reacting  sample. 

Figure  30.  Tridiagonal  matrixes  for  energy  equation  (Eq.25). 

Figure  31.  Screen  shoot  of  window  version. 

Figure  32.  Temperature  profile  in  a  non-reacting  cylindrical  sample  with  At  =  0.004. 

Figure  33.  Temperature  profile  in  the  cylindrical  sample  with  At  =  0.004. 

Figure  34.  Conversion  profile  in  the  cylindrical  sample  with  At  =  0.004. 

Figure  35.  Temperature  profile  in  a  reacting  cylindrical  sample  for  frequency  factor  of 
1.728E10S1  with  At  =  0.002. 

Figure  37.  Temperature  profile  in  the  reacting  cylindrical  sample  for  frequency  factor 
of  2.016el0s~I  with  At  =  0.002. 

Figure  38.  Residted  matrix  for  energy  equation. 

Figure  39.  Continuity  equation  matrix. 

Figure  40.  Dimensionless  pressure  profile  for  non-reacting  system  with  At  =  0.005. 

Figure  41.  Dimensionless  density  profile  for  non-reacting  system  with  At  =  0.005. 

Figure  42.  Dimensionless  temperature  profile  for  non-reacting  system  with  At  =  0.005. 

Figure  43.  Dimensionless  pressure  profile  with  At  =  0. 005. 

Figure  44.  Dimensionless  density  profile  for  50  bars  (725  psi)  with  At  =  0.005. 

Figure  45.  Dimensionless  temperature  profile  for  50  bar  (725  psi)  with  At  =  0.005. 
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Figure  46.  Dimensionless  velocity  profile  for  50  bar  (725  psi)  with  At  =  0.005. 

Figure  4  7.  Conversion  profile  for  50  bar  (725  psi)  with  At  =  0. 005. 

Figure  48.  Dimensionless  pressure  profile  with  At  =  0.004. 

Figure  49.  Dimensionless  temperature  profile  for  70  bars  (1000  psi)  with  At  =  0.004. 

Figure  50.  Dimensionless  density  for  70  bars  (1000  psi)  with  At  =  0.004. 

Figure  51.  Dimensionless  velocitv  profile  for  70  bars  (1000  psi)  pressure  with  At  = 

0.004. 

Figure  52.  Conversion  profile  for  pressure  70  bars  (1000  psi)  with  At  =  0.004. 

Figure  53.  Dimensionless  pressure  profile  with  At  =  0.003. 

Figure  54.  Dimensionless  temperature  profile  for  2000  psi  pressure  with  At  =  0.003. 

Figure  55.  Dimensionless  density  profile  for  2000  psi  pressure  with  At  =  0.003. 

Figure  56.  Dimensionless  velocity  profile  for  2000  psi  pressure  with  At  =  0.003. 

Figure  57.  Conversion  profile  for  2000  psi  pressure  with  At  =  0.003. 

Figure  58.  Dimensionless  velocity  profile  at  various  pressures. 
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Table  3.  Average  propagation  velocity  of  the  investigated  systems  under  unconfined 
conditions. 

Table  4.  Results  for  the  measurement  of propagation  velocity  of Al/CuO  as  a  function  of  weight 
percent  coating  on  the  aluminum  nano-powder.  The  coatings  used  were  Z6124  silane 
and  oleic  acid. 
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STATEMENT  OF  THE  PROBLEM  STUDIED 


Metal-based  energetic  materials,  under  certain  stimuli,  release  large  amounts  of  energy. 
They  are  essential  ingredients  in  rocket  propellants,  primers,  and  warhead  explosives.  However, 
rate  of  energy  release  for  these  metal-based  energetic  materials  are  rather  slow.  Metal-based 
energetic  materials  that  are  produced  in  nano-scale  are  potentially  attractive  due  to  their  rapid 
energy  release.  This  increased  performance  is  due  to  a  significant  increase  in  the  specific  surface 
area  of  all  involved  reactants.  As  the  specific  surface  area  increases,  the  number  of  contact  points 
between  the  reactants  also  increases.  Review  of  recent  literature  on  the  combustion  of 
pyrotechnic  materials  has  shown  that  there  is  a  significant  relation  between  reaction  rate  and 
particle  size.  Reducing  particle  size  results  in  increase  of  combustion  front  velocity  due  to 
reduced  diffusional  distances  between  reactant  particles.  Further,  reaction  rate  depends  on 
various  factors  including  particle  size  distribution  and  degree  of  intermixing.  Reaction  rates 
between  nanosize  aluminum  and  metal  oxides  can  be  significantly  greater  than  those  observed 
with  traditional  micron-size  thermite  powders.  One  class  of  nano-scale  metal-based  energetic 
materials  includes  systems  consisting  of  nanosize  aluminum  and  metal  oxide  as  an  oxidizer.  This 
type  of  reacting  system  is  known  as  metastable  intermolecular  composite  (MIC).  In  other  words, 
MIC  is  a  mixture  of  nano-scale  energetic  constituents  exhibiting  a  thermite  behavior.  Such 
reactions  occurring  between  metal  and  metal  oxide  powders  are  accompanied  by  the  generation 
of  high  temperatures  (>3000K).  Super-thermites,  formed  by  combining  aluminum  and  metal 
oxide  nanopowders  are  currently  incorporated  in  fonnulating  the  next  generation  of 
environmentally  friendly  primers  and  lead-free  matches.  They  may  also  find  an  application  in  air 
bag  initiators  and  inflators,  and  thermal  batteries.  Until  now,  general  combustion  characteristics 
of  MIC  reacting  systems  have  not  been  well  understood  therefore  this  research  study  was 
focused  on  better  understanding  of  combustion  front  propagation  characteristics  using  both 
experimental  and  theoretical  tools.  The  experimental  studies  were  focused  on  functionalization 
of  aluminum  nanopowder  in  order  to  reduce  the  aging  effect  of  moisture,  improve  its  dispersion 
in  organic  solvents,  and  improve  mixing  with  oxide  nanopowders.  The  effect  of  aluminum 
surface  modification  on  aging  and  reactivity  of  nanopowders  with  copper  oxide  were  studied 
extensively.  The  experimental  research  was  accompanied  by  numerical  studies  of  a  combustion 
front  propagation  in  the  presence  of  pressure  gradients  generated  during  the  reaction. 


SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

During  the  first  period  of  research  grant  October  1,  2002  -  December  31,  2003  the  following 
research  initiatives  have  been  undertaken: 

•  Functionalization  and  characterization  of  aluminum  nanopowders. 

•  Aging  studies  of  aluminum  nanopowders  at  different  relative  humidities. 

•  Investigation  of  reactivity  of  aluminum  with  copper  oxide  nanopowders  under 
unconfined  conditions. 

•  Effect  of  coating  of  aluminum  nanopowders  on  combustion  front  velocities. 

•  Mixing  studies  of  systems  consisting  of  binary  nanopowders. 
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•  Mathematical  model  describing  combustion  front  propagation  in  condensed  phase  system 
with  and  without  products  vaporization. 


1.  CHARACTERIZATION  OF  ALUMINUM  NANOPOWDERS 

Aluminum  is  very  reactive  metal  that  undergoes  quick  oxidation  and  hydration  at  ambient 
conditions.  Particularly,  the  nanosize  aluminum  particles  with  large  exposed  surface  can  react 
violently  with  oxygen  (nanosize,  pure  aluminum  metal  is  pyrophoric).  Preparation  of  oxide 
passivation  barrier  on  the  particle’s  surface  controls  this  process  and,  frequently,  the  nanosize 
aluminum  available  on  the  market  is  an  oxide-passivated  material. 

A  typical  TEM  image  of  passivated  A1  nanopowder  is  shown  in  Figure  1.  As  can  be  seen 
from  that  photograph,  a  thin  layer  (approximately  1.5-2. 5  nm)  of  aluminum  oxide  is  present  on 
the  aluminum  surface. 


Figure  1.  TEM  image  of  passivated  aluminum  nanopowder. 

Even  a  thin  layer  of  the  oxide  on  the  surface  amounts  for  large  reduction  of  the  reactive  A1 
content  of  the  passivated  material.  Quantitative  assessment  of  the  active  aluminum  in  the 
nanopowder  is  essential  for  designing  and  carrying  out  experiments  where  stoichiometry  of  the 
reactants  is  relevant.  Therefore,  our  research  studies  were  focused  on  quantitative  determination 
of  the  content  of  reactive  aluminum  in  passivated  A1  nanopowders.  Both,  TGA  and  volumetric 
method  were  used  in  these  studies. 

The  amount  of  reactive  aluminum  detennined  from  TGA  analyses  conducted  both  in  air 
and  oxygen,  is  apparently  lower  than  expected.  A  sudden  ignition  of  the  aluminum  sample  leads 
to  the  melting  of  the  aluminum  within  an  oxide  product.  This  melting  process  prevents 
completion  of  the  oxidation  process  by  entrapping  unreacted  aluminum  inside  the  macroscopic 
aluminum  oxide  droplet.  In  order  to  prevent  the  formation  of  such  droplets  in  next  experiments, 
submicron  alumina  powder  was  used  as  a  diluent.  It  was  determined  that  approximately  35  mg  of 
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AI2O3  was  needed  to  mix  with  5  mg  of  A1  nanopowder.  It  was  also  determined  that  an  additional 
10%  weight  increase  was  recorded  in  TGA  experiment  as  the  applied  temperature  increased  from 
800  to  1450°C.  In  order  to  ascertain  the  completeness  of  oxidation,  the  final  experiments  were 
carried  out  until  1450°C.  The  results  are  shown  in  Table  1. 

Table  1.  TGA  analyses  of  aluminum  nanopowders  obtained  from  NSWC/IH. 


A1  nanopowder 
(surface  area) 

Gas  Atmosphere 

Average  %  of  Reactive 
Aluminum  Content 
Temperature:  25-1450  °C 

1030001  -  (34.5  m2/g) 

Air 

68 

1-65-38  -(79.1  nr/g) 

Air 

48 

1-62-37  -  (50.3  m2/g) 

Air 

67 

Another  important  factor  affecting  TGA  analysis  of  A1  nanopowders  is  the  choice  of  a  reactive 
gas.  Normally,  air  is  used  in  most  combustion  experiments.  However,  it  was  determined  in  our 
experiments  these  research  studies  that  the  product  of  combustion  of  aluminum  in  air  is  not  pure 
aluminum  oxide  but  solid  solution  of  that  oxide  with  aluminum  nitride.  The  latter  component  is 
the  product  of  reaction  of  aluminum  with  nitrogen  in  air.  The  material  obtained  from  combustion 
in  air  is  alon  (trace  3  in  Figure  2). 


15  25  35  45  55  65  75 


two  theta  (degree) 


Figure  2.  X-ray  patterns  of passivated  aluminum  and  products  of  combustion  in  oxygen 

(alumina)  and  in  air  (alon). 

The  TGA  technique  is  not  very  suitable  for  measurement  of  reactive  aluminum  content  in 
polymer  coated  A1  nanopowders  due  to  simultaneous  oxidation  of  metal  and  polymer. 

Therefore,  another  method  based  on  wet  chemistry  (volumetric  method)  is  more  suitable  in  such 
situations.  The  volumetric  method  used  in  this  research  for  characterization  of  nanosized 
aluminum  powder  was  also  revised  to  introduce  improvements  in  accuracy  and  convenience.  The 
principle  concept  of  this  method  is  a  selective  reaction  between  aluminum  powder  and  water  in 
the  presence  of  sodium  hydroxide  (NaOH): 
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A1  +  3H20  ->  A1(0H)3  +  1.5H2T. 

By  measuring  volume  of  generated  hydrogen  one  can  infer  about  amount  of  equivalent 
aluminum  involved  in  the  reaction.  In  calculations  the  precise  determination  of  ambient  pressure 
and  temperature  is  required.  Saturation  with  water  of  nitrogen  gas  used  for  purging  the  apparatus 
is  also  needed  to  minimize  experimental  errors. 


Figure  3.  Schematic  of  the  glassware  for  the  determination  of  the  reactive  aluminum 
content  by  the  volumetric  method. 


Sample  cup 


Figure  4.  Schematic  of  the  stopper  used  in  neck  6.  Sample  cup  is  a  metal  piece  bonded 
to  the  end  of  the  glass  stopper. 

The  experimental  procedure  was  described  in  detail  in  the  previous  progress  report. 
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5.  AGING  STUDIES  OF  NANOSIZE  ALUMINUM  POWDERS. 


A  systematic  study  of  slow  oxidation  of  aluminum  nanopowders  in  presence  of  various 
partial  pressure  of  water  was  undertaken  to  develop  effective  measures  for  preventing  aging  of 
the  nanopowder  in  ambient  conditions. 

Sample  vials  containing  untreated  (as  received)  aluminum,  UFAL  1-65-38,  were  placed 
in  sealed  vessels  maintaining  relative  humidity  of  43,  75,  83,  and  97  percent.  The  sample  vials 
were  open  to  the  vessel  atmosphere.  The  vessels  were  placed  in  an  oven  maintained  at  40°C. 
See  Figure  5  below  for  a  schematic  of  the  humidity  vessel. 


Figure  5.  Schematic  of  an  aging  chamber. 

The  relative  humidity  (RFI)  inside  each  vessel  was  maintained  with  selected  saturated  aqueous 
salt  solutions  as  recommended  in  ASTM  Standard  E  104-85  Standard  Practice  for  Maintaining 
Constant  Relative  Humidity  by  Means  of  Aqueous  Solutions.  The  aqueous  salt  solutions  used 
were  potassium  carbonate,  sodium  chloride,  potassium  chloride,  and  potassium  sulfate.  See 
Table  2  below  for  the  salt  solutions  and  their  %RH. 

Table  2.  Salt  solutions  for  maintaining  humidity. 


k2co3 

NaCI 

KCI 

K2S04 

43%RH 

75%RH 

83%RH 

97%  RH 

Samples  of  the  aged  powders  were  tested  for  percent  reactive  A1  using  the  volumetric  method. 
Samples  were  taken  approximately  every  3  or  6  days  depending  on  the  humidity  the  sample  is 
exposed  to.  The  results  are  shown  in  Figure  6. 


9 


%  Reactive  Al  of  Powders 
Aged  at  40C  and  Specified  %RH 


Days  Aged 


Figure  6.  Aging  results  for  uncoated  Al  (UFA!  1-65-38). 

Untreated  aluminum  nanopowder  exposed  to  humid  air  undergoes  oxidation,  particularly  rapidly 
at  humid  conditions.  This  process  leads  to  drastic  reduction  of  the  reactive  aluminum  content.  In 
order  to  prevent  or  effectively  slow  down  deterioration  of  the  material,  application  of  various 
surface  coatings  were  considered. 

Aluminum  nanopowder  coating  experiments  are  currently  underway  utilizing  several  Dow 
Corning  silanes.  Two  samples  of  UFA1  1-65-38  aluminum  have  been  dispersed  with  3  wt% 
silane.  The  silanes  currently  being  examined  are  Dow  Coming  Z-6124  and  Z-6040.  The  coated 
Al  samples  were  aged  under  the  same  conditions  as  the  above  aging  experiment  to  determine 
their  effectiveness.  An  uncoated  aluminum  sample  was  aged  concurrently  with  the  coated 
samples  as  a  control.  The  coating  procedure  is  as  follows: 

Coating  Procedure: 

1 .  Al  sample  is  weighed  out  and  placed  in  a  vial. 

2.  The  required  weight  of  silane  is  added  to  the  vial.  The  silane  is  contained  in  a 
previously  prepared  stock  solution,  of  silane  and  methanol,  of  known  silane 
concentration.  A  calculated  weight  of  stock  solution  is  weighed  into  the  vial. 

3.  Methanol  is  used  to  dilute/disperse  the  Al  and  silane  in  the  vial. 

4.  The  vial  is  seal  and  sonicated  for  15  minutes. 

5.  After  sonication  the  dispersed  Al/silane/methanol  solution  is  poured  into  a  pan  and 
the  methanol  allowed  to  evaporate  leaving  behind  the  Al  and  silane. 

6.  The  Al/silane  powder  is  then  collected  from  the  pan. 

Results  for  the  aging  of  Al  coated  with  3  wt.%  Dow  Coming  Z-6124  and  Z-6040  silane  and 
uncoated  Al  are  shown  in  Figures  7-10. 
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43%  RH  Aging 

Coated  and  Uncoated  Powders  @  40C 


Figure  7.  Aging  results  for  coated  and  uncoated  Al  at  43  %RFl. 


75%  RH  Aging 

Coated  and  Uncoated  Powders  @  40C 


Days  Aged 


Figure  8.  Aging  results  for  coated  and  uncoated  Al  at  75  %RH. 
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84%  RH  Aging 

Coated  and  Uncoated  Powders  @  40C 


Figure  9.  Aging  results  for  coated  and  uncoated  Al  at  84  %RH. 


97%  RH  Aging 

Coated  and  Uncoated  Powders  @  40C 


Days  Aged 


Figure  10.  Aging  results  for  coated  and  uncoated  Al  at  97  %RH. 


Silane  coating  of  aluminum  nanopowders  results  in  several-fold  decrease  of  aging  rate  as 
compared  to  an  untreated  material.  Both  silanes  are  similarly  effective  at  RH  <85%.  The  Z6124 
silane  appeared  as  more  effective  coating  than  Z6040  silane  only  at  very  high  relative  humidity. 

Work  is  underway  to  determine  the  amount  of  silane  (or  other  coating  agents)  that  has  bonded 
(chemically  or  via  adsorption)  to  the  surface  of  the  nanopowders.  The  experiments  will  make 
use  of  the  coating  agent  solution  concentration  before  and  after  the  powders  have  been  dispersed 
in  the  solution.  To  eliminate  the  coating  agent  that  maybe  adsorbed,  the  powders  will  then  be 
washed  in  pure  solvent  and  the  solvent  analyzed  to  determine  if  the  coating  agent  is  present. 
Currently  Dow  Corning  Z-6124  Silane  is  being  examined  by  UV-Vis.  UV-Vis  is  being  used  for 
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it’s  ease  and  short  analysis  time.  The  unique  structure  of  this  silane  allows  for  analysis  is  the  UV 
range.  Z-6124  is  phenyltrimethoxysilane.  The  phenyl  group  shows  relatively  strong  absorption 
in  the  260-275  nm  range. 

Although  the  43  %  RH  samples  react  slower  than  those  at  higher  humidities,  a  nearly  60 
percent  decrease  in  the  reactive  aluminum  content  in  50  days  is  still  significant.  A  relative 
humidity  of  43  percent  is  not  an  unusual  ambient  reading.  If  not  adequately  sealed  from  the 
environment,  a  dramatic  decrease  in  the  reactivity  of  aluminum  nano-powders  or  reactive 
systems  utilizing  aluminum  will  be  experienced. 

It  was  found  that  the  curves  in  above  presented  figure  initially  show  a  flat  slope  with  little 
decrease  in  the  percent  reactive  aluminum.  A  steep  slope  corresponding  to  a  dramatic  decrease 
in  the  percent  reactive  aluminum  follows  the  flat  portion  of  the  curve.  In  the  flat  region  of  the 
curve  it  is  believed  that  hydration  of  the  existing  oxide  layer  is  taking  place.  Once  the  existing 
layer  has  been  hydrated,  reaction  of  the  core  (unreacted  aluminum)  begins,  which  is  represented 
by  the  sharp  increase  in  the  slope  of  the  curves. 

Similar  results  for  hydration  were  shown  by  Grimblot  and  Eldridge  [1],  which  they 
tenned  an  incubation  period  (IP).  In  their  research  the  oxidation  of  A1  was  monitored  by  weight 
change  of  aluminum  samples  immersed  in  aqueous  solutions.  Aluminum  oxide  layers  were 
formed  on  the  aluminum  samples  prior  to  immersion  in  solution.  It  was  found  that  when  placed 
in  aqueous  solution,  aluminum  samples  with  oxide  layers  formed  in  a  dry  environment  show  a 
period  of  no  weight  change,  or  an  IP,  followed  by  a  fast  rate  of  weight  change.  However, 
samples  prepared  with  a  hydrated  oxide  layer  did  not  exhibit  the  IP  but  immediately  began  to 
show  a  fast  rate  of  weight  change.  Grimblot  and  Eldridge  believe  that  the  hydrated  AI2O3 
structure  is  more  open  allowing  for  easier  movement  of  diffusing  species  (i.e.  water). 

Similar  results  were  found  when  aging  samples  of  50  nm  Nanotechnologies  and  51  nm 
Technanogy  aluminum  in  97  %RH  at  40  °C.  The  results  are  shown  in  Figure  11  below.  The 
Nanotechnologies  sample  shows  a  period  of  little  to  no  change  in  percent  reactive  aluminum 
until  around  26  hours,  while  the  Technanogy  sample  shows  immediate  decreases  in  percent 
reactive  aluminum  after  8  hours.  Samples  of  both  powders  were  analyzed  for  water  content 
using  TGA.  It  was  found  that  the  Technanogy  sample  contained  between  7  and  8  wt%  water. 
The  Nanotechnologies  sample  only  contained  around  3  wt%  water.  It  is  likely  that  the  oxide 
layer  of  the  Technanogy  sample  is  more  hydrated  compared  to  the  Nanotechnologies  sample. 
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Time  (hours) 


Figure  11:  Results  for  aging  of  Nanotechnologies  and  Technanogy  aluminum  samples  in 
97  %RH  at  40  °C.  Lines  shown  are  trendlines  only  and  each  data  point  is 
an  average  of  two  volumetric  method  measurements. 


5.  MIXING  OF  NANOSIZE  ALUMINUM  AND  TITANIUM  DIOXIDE  POWDERS. 

The  objective  of  the  experiment  was  twofold: 

•  To  study  the  formation  of  agglomeration  and  the  quality  of  mixing; 

•  to  compare  hexane  and  ethanol  as  solvent  for  nanosized  powder  systems. 

Aluminum  and  titanium  dioxide  powders  were  chosen  as  the  model  reacting  system. 
Assuming  the  same  particle  size  of  aluminum  and  titanium  dioxide  powders  and  considering  the 
difference  in  the  density  of  these  two  materials  (A1  2.7  g/cm3  and  TiC>2  4.26  g/cm3)  a  mixture  of 
such  two  powders  in  weight  ratio  equal  to  1 : 1  molar  ratio,  will  has  the  property  to  occupy  the 
similar  volume  by  each  of  the  components.  Volume  calculation  based  on  the  stated  assumptions 
result  in  a  Aluminum:  Titanium  dioxide  ratio  of  41:59.  This  weight  ration  was  kept  constant  for 
all  experiments. 

Preparation  of  the  samples. 

Absolute  ethanol  and  n-hexane  from  stock  was  used  as  a  solvent.  Two  weight  %  of 
Sodium  dioctyl  sulfosuccinate  (“Aerosol  OT  100%  Surfactant”  from  Cytec  Industries  Inc.)  was 
added  as  surfactant  to  each  solution. 

0.00269g  of  A1  and  0.00599g  of  TiCb  was  added  each  solution.  The  dispersions  were  mixing  in 
an  ultrasonic  bath  for  20  minutes.  After  mixing,  powder  mixtures  were  dried  on  aluminum  plate 
and  sample  were  taken  with  an  adhesive  carbon  tape.  The  carbon  tapes  were  placed  on  specimen 
holders  and  carbon  coated. 
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SEM/Microprobe  analyses. 

SEM  and  Microprobe  analyses  followed:  First,  secondary  electron  images  (which  have  a  good 
resolution)  and  backscattered  electron  images  (which  provide  better  contrast)  were  taken  at 
10000X.  (Figures  12,13  and  14,15).  However,  those  images  do  not  take  the  surface  roughness 
into  account. 


Figure  12.  Al-TiO 2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (SE-image  10000  X). 


Figure  13.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (B SE-image  10000  X). 


Figure  14.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (SE-image  10000  X). 
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Figure  15.  Al-TiO 2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (BSE-image  10000  X). 

After  that,  so  called  cameo  images  were  taken  as  first  screening.  The  Cameo  Image  technique  is 
a  half-quantitative  imaging  technique  developed  by  Oxford  Inc.,  combining  SE-image  signal 
with  energy  dispersive  x-ray  analysis  (see  Figures  16  and  17).  However  just  the  overall  energy 
values  of  a  sample  are  taken  into  account  and  not  the  specific  peak  of  element.  In  our  case,  the 
energy  value  of  aluminum  is  lower  than  the  energy  peaks  of  titanium,  so  we  can  paint  the  dots, 
which  have  lower  energy  value  (Al)  red  and  dots,  which  exhibit  higher  energy  value  blue  (Ti). 
The  green  color  represents  mixed  areas. 


Figure  16:  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (Cameo-image  10000  X). 
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Figure  17.  Al-TiO 2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (Cameo-image  10000  X). 


The  microprobe  elemental  analysis  were  conducted  as  mapping  tests  and  as  line  scans,  by  using 
the  strongest  energy  disperse  peaks  of  aluminum  and  titanium. 

The  mapping  was  taken  at  a  magnification  of  10000X.  Figures  18  and  19  show  SE-images,  BSE- 
images,  elemental  A1  -images,  and  elemental  Ti  -images.  This  overview  allows  distinguishing 
between  surface  roughness  (shadow  effects)  and  agglomerations. 


Figure  18.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (Elemental  Mapping  10000  X SE-BSE-Al  and  Ti  mapping  image). 
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Figure  19.  Al-TiO 2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  ( Elemental  Mapping  1 0000  X SE-BSE-Al  and  Ti  mapping  image). 

The  line  scan  pick  up  the  presence  of  both  elements  at  the  same  time  and  allows  for  a 
quantitative  comparison  of  the  x-ray  energy  signals,  which  can  be  relates  to  the  weight  and 
atomic  percentage  of  the  elements  present  in  the  sample.  Examples  of  the  line-scan  are  shown  in 
Figures  20  and  2 1 . 


Figure  20.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  (Three  line  scans  of  10  pm  at  10000  X). 


18 


Figure  21.  Al-TiO 2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant  ( Three  line  scans  of  10  pm  at  10000  X). 


It  was  interesting  to  notice  that  long  in  time  line  scans  (100s)  lead  to  reaction  of  the  samples  (see 
Figures  22  and  23),  whereas  the  mapping  or  the  cameo  analyses  (scanning  a  whole  area)  leaves 
the  sample  as  it  is.  It  should  be  further  noted  that  that  the  sample  prepared  with  abs.  ethanol 
(Figure  23)  shows  a  bigger  effect  than  the  sample  prepared  with  n-hexane  (Figure  24). 

Obviously  the  intensity  (energy  /area*time)  of  the  electron  beam  was  sufficient  to  promote  a 
reaction.  The  subsequent  spot  analysis  conducted  at  the  sample  which  was  prepared  with  ethanol, 
proves  that  there  are  differences  in  the  composition  of  the  areas  where  lines  scan  were  made 
(Figure  23:  spot  #1  and  #4)  and  where  no  lines  scans  were  made  (Figure  23:  spot  #2  and  #3). 
However,  how  these  changes  finally  affect  the  results  was  not  further  investigated  due  to  safety 
reasons. 


Figure  22.  Al-TiO  2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctvl  sulfosuccinate  as 
surfactant.  Sample  after  three  line  scans  of  10  pm  at  10000  X. 


19 


Figure  23.  Al-TiO 2-mixture  prepared  in  absolute  ethanol  with  sodium  dioctyl  sulfosuccinate  as 
surfactant.  Sample  after  three  line  scans  of  10  pm  at  10000  X.  The  X-ray  spot 
analyses  at  four  points  for  aluminum  and  titanium  show  the  following  results  #1  (Al: 
44.7  Ti:  55.3),  #2  (Al:  48.4  Ti:  51.6),  #3  (Al:  48.8  Ti:  51.2),  #4  (Al:  40.8  Ti:59.2). 


Figure  24.  Al-TiO  2-mixture  prepared  in  n-hexane  with  sodium  dioctyl  sulfosuccinate  as 
surfactant.  Sample  after  the  three  line-scans  of  10  pm  at  10000  X. 

The  mixing  technique  with  ethanol  seems  to  work  better  based  on  the  qualitative  results. 
However,  the  resulting  mixture  seems  to  be  also  more  reactive.  The  next  step  of  these  mixing 
experiments  will  be  the  testing  of  combustion  velocity  of  those  powder  mixtures. 


5.  MEASUREMENTS  OF  PROPAGATION  VELOCITIES. 

Combustion  front  propagation  velocity  for  unconfined  reactive  powders  is  one  of  the 
indirect  measures  of  the  inhomogeneous  reaction  rate  and  the  reactant  affinity.  Propagation 
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velocities  in  Al-CuO  and  AI-M0O3  reacting  systems  were  found  to  have,  very  high  combustion 
front  velocities.  Moreover,  due  to  a  rapid  expansion  of  gases  generated  at  very  high  temperatures 
in  Al-CuO  and  AI-M0O3  systems,  unreacted  powders  ahead  of  the  combustion  front  are 
frequently  expelled  from  the  bum  tray’s  cavity.  Because  of  this  displacement,  the  measured 
combustion  velocities  might  be  erratic.  In  order  to  countenneasure  this  effect  series  of  perforated 
baffles  was  inserted  into  the  trough  to  minimize  powder  displacement.  The  propagation 
velocities  in  the  Al-CuO  system,  measured  without  the  use  of  baffles,  varied  from  100  m/s  to 
1000  m/s.  However,  when  baffles  were  used,  the  combustion  front  velocities  in  the  same  system 
were  consistent  and  significantly  lower  (250  +  25  m/s).  Table  3  shows  the  average  combustion  front 
velocities  in  Al-CuO  reacting  systems  in  an  air  atmosphere. 

Table  3.  Average  propagation  velocity  of  the  investigated  systems  under  unconfined  conditions. 


Reacting  System 
(Nanosize  Powders) 

Atmosphere 

Combustion  Front 
Velocities  [m/s] 

A1  50  nm(Nanotechnologies  Inc) 
and  CuO  (Nanophase  Tech) 

Air 

416  +  25 

A1  80  nm(Nanotechnologies  Inc) 
and  CuO  (Nanophase  Tech) 

Air 

500  ±25 

Direction  of  propagation  of  the  combustion  wave  front 


t=0.1  fJS 


t=  50  fj  s 


t=  100  //  s 


t=  150  fj  s 


t=  1 60  n  s 


t=  170  fj  s 


t=  200  n  s 


Figure  25.  Images  of  burning  loose  nanopowders  of  AI-M0O3  system  recorded,  with  perforated 

baffles,  using  high-speed  camera. 
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Figure  25  shows  typical  images  of  a  combustion  front  propagation  in  the  AI-M0O3 
system.  The  recording  rate  was  8,000  fps  with  a  shutter  speed  of  1/128,000  s.  Bright  plumes, 
shown  in  the  images,  are  likely  composed  of  metal  vapor  and  alumina  particulates  when  the 
combustion  was  carried  out  in  an  inert  atmosphere.  When  the  thermite  reaction  took  place  in  air, 
a  significant  plume  was  formed  due  to  the  secondary  reaction  between  generated  vapor  and 
gaseous  oxygen. 

Silane  and  oleic  acid  coatings  prove  to  be  effective  at  preventing  the  aging  of  aluminum 
nano-powders  but  it  was  not  known  what  affect  they  would  have  on  the  reactivity  of  MIC 
reactant  mixtures  (specifically  aluminum/copper  oxide)  containing  coated  aluminum.  Samples 
of  aluminum  containing  0,  1,  3,  5,  10,  and  15  wt%  silane  were  mixed  with  CuO  and  tested  for 
propagation  velocity  in  open  burn  tray  measurements.  The  results  can  be  seen  in  Table  4  and 
Figure  26  below. 


Figure  26.  Results  for  the  measurement  of propagation  velocity  of Al/CuO  as  a 
function  of  weight  percent  coating  on  the  aluminum  nano-powder. 
The  coatings  used  were  Z6124  silane  and  oleic  acid.  Each  data 
point  is  the  average  of  4  measurements. 
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Table  4.  Results  for  the  measurement  of propagation  velocity  of Al/CuO  as  a 

function  of  weight  percent  coating  on  the  aluminum  nano-powder. 
The  coatings  used  were  Z6124  silane  and  oleic  acid. 


Z6124  Silane 

oleic  acid 

Weight  % 

*Propagation 

Std  dev 

*Propagation 

Std  dev 

Coating 

Velocity  (m/s) 

(m/s) 

Velocity  (m/s) 

(m/s) 

0 

528 

27 

528 

27 

1 

538 

30 

589 

49 

3 

500 

46 

673 

67 

5 

576 

47 

747 

56 

10 

580 

70 

780 

14 

15 

*Each  data  point  is  an  average  of  4  measurements. 


It  is  shown  as  oleic  acid  concentration  is  increased  up  to  10  weight  percent  there  is  a 
dramatic  increase  in  the  propagation  velocity.  In  all,  there  is  almost  a  50  %  increase  in 
propagation  velocity  for  the  10  weight  percent  sample  compared  to  the  uncoated  sample.  Data 
appears  to  be  leveling  off  between  5  and  10  weight  percent  and  at  15  weight  percent  oleic  acid 
ignition  does  not  take  place.  This  indicated  there  is  an  optimal  concentration  around  10  weight 
percent. 

The  Z6124  silane  coated  samples  show  a  similar  but  less  dramatic  trend.  In  fact  if  the 
standard  deviation  is  taken  into  account  the  velocities  measured  cannot  be  considered  different 
with  any  reasonable  confidence.  As  was  seen  with  oleic  acid,  the  15  weight  percent  sample  also 
was  unable  to  ignite. 

Propagation  velocity  is  increased  due  to  enhanced  mixing  of  the  reactants  with  increased 
coating  concentration.  It  is  now  known  for  50  mn  aluminum  that  Z6124  silane  in  excess  of  3 
wt%  is  not  chemically  bonded  to  the  aluminum  and  would  disperse  in  the  mixing  solvent 
(hexane).  The  excess  silane  then  acts  as  a  surfactant  improving  mixing  between  both  reactants. 
The  bonded  silane  acts  as  a  surfactant  as  well,  helping  to  suspend  the  aluminum  particles  in 
solution.  However,  as  the  coating  concentration  increases  the  coating  layer  begins  to  thicken. 
The  coating  creates  an  additional  barrier  between  reactants  that  must  be  removed  for  reactants  to 
come  into  contact.  This  increase  in  the  barrier  works  in  opposition  to  the  enhancement  of 
mixing,  and  at  15  wt%  coatings  are  in  too  large  a  concentration  and  prevent  contact  between  the 
aluminum  and  copper  oxide  reactants.  If  the  reactants  are  not  in  contact  ignition  is  not  possible. 

Furthermore  the  oleic  acid  shows  improvement  over  Z6124  silane  for  three  reasons.  First 
physisorption  is  more  likely  to  take  place  for  oleic  acid  than  chemisorption.  The  physisorb  oleic 
acid  molecules  disperse  in  the  mixing  solvent  and  can  contact  both  the  aluminum  and  copper 
oxide  particles  equally.  The  second  reason  is  the  long  18-member  oleic  acid  tail  has  greater 
molecular  interaction  with  the  straight  chain  hexane  molecule  than  the  phenyl  group  on  the 
Z6124  Silane.  The  enhanced  interaction  makes  oleic  acid  a  better  surfactant  improving  mixing 
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of  the  reactants.  The  final  reason  is  that  after  hydrolysis  of  the  silane  to  form  silanol,  self¬ 
condensation  of  silanol  can  occur  to  form  siloxane  linkages.  The  formation  of  the  siloxane 
linkages  may  create  an  additional  barrier  (in  the  form  of  a  siloxane  layer)  that  needs  to  be  broken 
down  for  reaction  to  occur. 

To  help  support  the  above  hypothesis  of  improved  mixing,  a  second  set  of  experiments 
were  performed  using  oleic  acid  and  1-6341  silane,  designated  n-octyltriethoxysilane.  1-6341 
has  a  straight  8-member  hydrocarbon  chain  as  its  organic  group.  It  was  assumed  that  1-6341 
would  have  improved  mixing  capability  over  Z6124  because  of  the  long  carbon  chain,  which 
should  have  stronger  interaction  with  hexane.  In  this  experiment  the  aluminum  powders  were 
not  coated  prior  to  mixing.  The  “coating”  was  added,  at  various  concentrations,  directly  to 
hexane  during  the  mixing  of  the  Al/CuO  reactants.  The  amounts  of  coating  added  were  exactly 
the  same  as  the  amounts  present  in  the  previous  experiments.  The  mixed  powders  were  dried 
and  the  coating  remained  on  the  powders.  The  results  were  as  expected,  as  coating  concentration 
increase  propagation  velocity  increased  for  both  oleic  acid  and  1-6341  silane.  Comparing  these 
results  to  the  previous  experiments,  the  percent  increase  in  propagation  velocity  from  lowest 
concentration  to  highest  concentration  (10  wt%)  is  still  close  to  50  %  for  oleic  acid.  Also  as 
expected,  1-6341  silane  was  an  improvement  over  Z6124  silane.  Comparing  the  results  of  the  0 
and  10  wt%  concentrations,  1-6341  silane  had  a  40  %  increase  in  propagation  velocity  compared 
to  less  than  a  10  %  increase  for  Z6124  silane.  The  results  are  shown  in  Figure  27  below. 
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Figure  27.  Propagation  velocity  as  a  function  of  coating  concentration. 

Coatings  were  added  at  the  time  of  mixing  of  AI/CuO  reactants. 

Effect  of  Coating  on  Ignition  delay 

It  was  initially  expected  that  coatings  on  aluminum  nano-powder  would  add  an  additional 
barrier  that  needed  to  be  removed  for  reactants  to  contact  and  reaction  to  progress.  Experiments 
were  performed  in  which  the  weight  percent  coating  on  aluminum  was  varied  and  the  ignition 
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delay  time  of  Al/CuO,  mixed  using  these  coated  powders,  was  measured.  Both  oleic  acid  and 
Z6124  silane  were  tested. 

It  is  observed  that  the  silane  coating  affects  the  ignition  delay  to  a  greater  degree  than 
oleic  acid.  The  ignition  delay  was  2.5  times  greater  for  the  7  wt%  coating  of  silane  on  aluminum 
than  the  uncoated  aluminum  sample.  In  comparison,  ignition  delay  of  the  10  wt%  oleic  acid 
coated  aluminum  was  only  1 .4  times  greater  than  uncoated  aluminum.  The  results  are  shown  in 
Figure  28  below.  The  overall  results  were  as  expected  however;  coatings  on  aluminum  nano¬ 
powder  increase  the  ignition  delay  of  Al/CuO.  The  coatings  add  an  additional  barrier  that  must 
be  removed  for  reaction  to  initiate.  Just  as  the  aluminum  oxide  layer  must  be  removed  for 
reactants  to  come  into  contact,  the  coating  must  be  removed.  Additional  energy  is  required  to 
remove  the  coating  layer. 


Figure  28.  Al/CuO  ignition  delay  results  for  coated  aluminum  sample. 
Coatings  studied  were  Z6124  Silane  and  oleic  acid. 


5.  MATHEMATICAL  MODELING  OF  CONDENSED-PHASE  COMBUSTION. 

Two  different  models  for  analyses  of  combustion-wave  propagation  in  condensed  phase 
reactions  were  examined  in  these  research  studies: 

I.  Condensed-phase  reactions  without  significant  vaporization  of  reactants  and  products; 

II.  Condensed-phase  reactions  with  significant  vaporization  of  reactants  and  products. 
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The  objective  was  to  investigate  on  1-D  model  for  both  condensed  phase  reactions  with  and 
without  significant  vaporization. 

Mathematical  model  for  condensed-phase  reaction  without  significant  vaporization  of 

reactants  and  products 

In  the  first  model,  the  mathematical  model  was  done  in  a  cylindrical  reacting  sample  along  which 
the  self-propagating  high  temperature  combustion  wave  occurs.  Figure  29  shows  the  cylindrical 
sample  over  which  the  material  and  energy  balance  were  modeled  [2], 


Figure  29.  Schematic  of  the  cylindrical  reacting  sample. 


The  material  balance  over  the  cylindrical  reacting  sample  mixture  is  given  by: 

=  (p{ri’T)  (1) 

at 


The  energy  balance  over  the  reacting  sample  mixture  is: 


k.dF 

dz 


r  dT  v 
PsCps  ~z~  = - — 

ot  dz 


+  ■ 


PSW , 


sr  frac 

M 


(-a ujMnJ)) 


(2) 


Dimensional  analysis  for  condensed-phase  reaction  without  significant  vaporization  of 

reactants  and  products 

Due  to  computational  advantages  it  was  decided  that  Equations  1  and  2  should  be 
rewritten  into  a  dimensionless  form.  The  heat  release  function  (p{rj,T )  is  typically  represented 
as: 

=  (3) 

where: 

f{r])  =  \-rjp  (4) 
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Eq.5  could  be  approximated  using  Frank-Kamenetskii  approximation  [3]  written  as: 


(  E\  (  E)  f E(T-T ']) 

exp - =  exp - exp  — - — - — 

l  RT)  \  RTe)  \  RTC2 

in  which  Tc  represent  the  combustion  or  ignition  temperature.  It  was  decided  that  the 
dimensionless  temperature  had  to  be  driven  from  Eq.6  and  could  be  defined  as: 

„  T-T 


Substituting  Eq.7  into  Eq.6  gives: 


E  \  E  E 

exp - =  exp - exp  - 0 

l  RT  RT  RT 


It  was  decided  that 


tH  =  exp 


Thus,  Eq.8  become: 


exp - =  9T 

RT 


Solving  for  T  from  Eq.7  gives: 

T=Tc(e+i) 

Furthermore,  the  dimensionless  length  was  defined  as: 


Thus,  the  differential  operators  for  Eq.l  1  and  12  are: 


clT  =  Td 0 


dz  =  Ld% 


(13,  14) 


Rewriting  Eq.2  in  term  of  dimensionless  temperature  and  length  result  into: 


PsCpsTc  —  = 


d 0  kTr.  82e  PsWfrackQ 


dt  L2  d%2 


{- ah  PAn 


Simplifying  Eq.l 5  gives: 
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(16) 


a© 


k,  a2©  wfracPsko 


+■ 


*  P.Cp.L  9?  P,Cp,MT ; 


The  dimensionless  time  and  its  differential  operator  were  defined  as: 


t 

r  =  — 

t 


dt  =  t  dr 


(17,  18) 


where: 


t  = 


PsCPsL2 


Eq.  16  could  be  rewritten  as: 


(19) 


dr  d%2  V  pA  ’ 
where: 


¥  = 


k„MT ' 


(20) 


(21) 


The  dimensionless  fonn  for  mass  balance  is: 


(22) 


Boundary  and  initial  conditions 

Initial  conditions  are: 

z-  =  0  0<£<1:  0  =  T"  ~T(:  ,  7  =  0 

The  boundary  conditions  can  be  represented  as  follows: 


r  >  0  £  =  0:  0  =  0 

|  =  1:  —  =  0 

d% 


Crank-Nicolson  Implicit  Scheme  Method 

It  was  decided  to  use  a  numerical  discrete  Crank-Nicolson  implicit  scheme  for 
approximation  of  time  and  spatial  derivatives  in  both  energy  and  mass  balances  equations  [4-8]: 
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(23) 


Furthermore,  Eq.20  after  discritization  becomes: 


©;,+1  -0; 

Ar 


0 "V  -  20 "+1  +  0* 


n+ 1 


1+ 1 


\«+l 
2- 1 


2(A^)2 
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-2@"  +©m 


2(Af)2 


+  +, 


(24) 


Matrix-free  Approximation 

In  order  to  develop  a  code  for  the  generalized  equations,  it  was  decided  that  a  matrix 
form  for  Eq.24  is  necessary  to  solve  the  above  system  of  equations.  It  was  further  decided  that 
the  conversion  (Eq.23)  has  to  be  calculated  at  later  time  steps  to  avoid  complication  of  Eq.24. 
Thus,  Eq.23  was  written  explicitly.  Simplifying  Eq.24  gives: 


.ref  +B"Q 


72+1 


+  C70;;/  =£>," 


(25) 


Where: 

a:  =  - 


Bn  =1  + 
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At 
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2Ag 
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®w  -2©''  +©”-! 


2A£2 


+  A  AY, 


Applying  grids  for  i  -  1,2 to  Eq.25  result  into  a  matrix  (see  Fig. 30). 
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Figure  30.  Tridiagonal  matrixes  for  energy  equation  (Eq.25). 
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This  is  a  tridiagonal  system  i.e.  a  system  of  equations  with  finite  coefficient  only  on  the  main 
diagonal  (the  Bi ' s  ),  the  lower  diagonal  (the  /f ' s  ),  and  the  upper  diagonal  (the  Cl' s). 

A  standard  method  for  solving  a  system  of  linear,  algebraic  equations  is  Gaussian  elimination. 
Thomas’  algorithm  is  essentially  the  result  of  applying  Gaussian  elimination  to  the  tridiagonal 
system  of  equations.  A  code  was  written  in  FORTRAN  language  to  solve  Eq.23  and  the  above 
matrix.  The  developed  code  was  then  incorporated  into  a  window  format  using  Visual  basic  (See 
Fig. 31).  The  code  for  the  above  equations  was  named  1-D  implicit  (open  system)  as  seen  on 
Fig.3 1 .  However,  more  testing  is  necessary  on  the  window  version  to  make  it  user  friendly. 


Figure  31.  Screen  shoot  of  window  version. 


Numerical  Results 

Initially,  the  model  was  tested  for  non-reacting  cylindrical  sample.  The  frequency  factor, 
length  of  non-reacting  cylindrical  sample,  integral  time,  ignition,  and  initial  temperatures  were 
set  to  Os'1,  0.0254m,  0.004,  2000K,  and  300K  respectively.  Figure  32  shows  the  temperature 
profile  along  the  length  of  the  non-reacting  sample  for  various  times.  It  can  be  seen  that  in  non¬ 
reacting  cylindrical  sample,  the  sample  was  only  being  heated. 
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Figure  32.  Temperature  profile  in  a  non-reacting  cylindrical  sample  with  Ar  =  0.004. 

Secondly,  a  reacting  cylindrical  sample  was  considered  in  which  the  frequency  factor  was  set  to 
1.44E10S'1  with  the  rest  of  the  parameters  being  the  same  as  in  non-reacting  sample.  The  code 
generates  the  temperature  and  conversion  profile  along  the  cylindrical  sample  (see  Figure  33). 


Figure  33.  Temperature  profile  in  the  cylindrical  sample  with  Ar  =  0.004. 
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Figure  34.  Conversion  profile  in  the  cylindrical  sample  with  At  =  0.004. 


To  further  the  investigation,  the  length  of  the  cylindrical  sample  was  kept  at  0.0254  m  and  the 
pre-exponential  factor  was  increased  to  1.728E10S'1  and  2.016E10S'1  respectively.  These  changes 
were  made  to  investigate  the  effect  of  pre-exponential  factor  along  the  length  of  cylindrical 
sample.  The  dimensionless  temperature  profiles  that  show  these  changes  are  shown  in  Figures  35 
and  36. 
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Figure  35.  Temperature  profile  in  a  reacting  cylindrical  sample  for  frequency  factor  of 
1.728E10S1  with  Ar  =  0.002. 

It  can  be  seen  from  Figures  35  and  36  that  as  the  pre-exponential  factor  is  increased  the 
propagation  become  faster.  Figure  37  shows  the  temperature  profile  for  pre-exponential  factor  of 
2.016el0s'1. 


Figure  37.  Temperature  profile  in  the  reacting  cylindrical  sample  for  frequency  factor 
of2.016el0s~J  with  Ar  =  0.002. 
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Mathematical  model  of  combustion  front  propagation  under  the  presence  of  uniaxial  gas 
pressure  gradient 

In  many  cases  temperatures  generated  during  the  combustion  are  very  high  causing 
partial  or  complete  vaporization  of  product(s).  A  mathematical  model  describing  propagation 
characteristics  in  a  condensed-phase  system  with  significant  pressure  generation  behind  a 
combustion  front  was  developed  and  explored. 

A  one  dimensional  adiabatic  model  is  presented  below.  The  subscript  s  denotes  the  solid  phase 
and  the  subscript  g  denotes  the  gas  phase. 


dT 


dT 


pCp~ 7  +  vgPgCpg  —  =  — 
dt  dz  dz 


dr 

k,  — 
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dz 


+  - 


PSW, 


frac 


M 


(30) 


where: 

pCp  =  (l  -  <X>)psCps  +  OpgCpg 
Cpg  =  a  +  f3T  +  yT2 

The  material  balance  over  cylindrical  sample  is  given  by: 


Where  (p{i j  ,  T )  was  given  by  Eq.3  rewritten  as: 

(pivJ)  =  f(vp)f{T) 


The  continuity  equation  over  the  cylindrical  sample  was  defined  as  [9]: 

dPg  ,  kvJ_Q 

dt  dz 


(31) 

(32) 


(34) 


(3) 


(35) 


It  was  decided  that  the  pressure  drop  along  the  cylindrical  sample  could  be  evaluated  using 
Ergun  equation  given  by  [10]: 


dP 

dz 


G 

r(i -<*>)") 

"150(1-®)//  ' 

PgScDp 

l  J 

i  1  .  /  JL/ 

Dp 

(36) 


where: 


G  =  pg®vg  (37) 

Substituting  Eq.3 7  into  Eq.3 6  gives: 
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dz  gcD 


In  Eq.38  the  only  parameter  that  changes  with  pressure  is  the  density  of  the  gas.  It  was  further 
decided  that  the  mass  flow  rate  at  any  point  down  the  reactor  is  constant.  Furthermore,  it  was 
decided  that  the  gas  used  in  this  system  is  an  inert.  Thus,  the  density  of  the  gas  was  defined  as: 

CtYp} 


Pg  Pg0[  T  lip 


Substituting  Eq.39  for  density  into  Eq.38  gives: 


dz  gcD 


(l-®)f  150(i-®V  +  j  ,75  feY  L 


-2 -  —  ®v 

8o{  tap)  s 


In  addition  to  Eq.39  that  relates  the  density,  pressure,  and  temperature,  it  was  further  decided  that 
the  gas  used  in  the  system  behave  ideally.  Further  relation  between  density,  pressure,  and 
temperature  was  given  by  ideal  gas  law  written  as: 


In  addition  to  the  above  equations,  it  was  decided  that  the  porosity  could  be  formulated  as 

®=®c+A®?7/;  (42) 

Where  ®0  and  A®  are  the  initial  and  change  of  porosity  respectively. 

Dimensional  analysis  of  the  condensed-phase  reaction  with  significant  vaporization  of 
reactants  and  products 

Defining  the  dimensionless  variables,  the  heat  released  function  given  by  Eq.3  had  been 
previously  defined  in  a  dimensionless  form: 

•p(trT)  =  *./(>?,  )»(<M)  («) 

Thus,  the  temperature,  length,  and  their  differential  operators  are  given  by 

T  =  Tc  (0  +  1)  £  =  j-  dT  =  Tcd&  dz  =  Ld%  (44) 


dT  =  TdQ 


dz  =  Ldd, 


Rewriting  the  energy  equation,  Eq.30,  into  dimensionless  length  and  temperature  results: 


PCpTc^~  + 
ot 


50  vsPgCpgTc  50  Tc  5  (,  50  )  PsWfrac 
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L  L2  d%) 


Simplifying  Eq.45  gives: 
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pCp  =  (1  -  Oa)psCps  +  ®0pgCpg 

The  dimensionless  velocity,  density  and  time  were  defined  as: 

Ps 


t  = 


pCpL 
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T  =  — 
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P  =' 


Pg 


Thus,  Eq.46  could  be  rewritten  in  dimensionless  form  as: 

S@+A0S  =  S^  ()„(<,-,) 
dz  dd,  dZ;  V; 

where: 


Y  = 


M-kT 


dt  =  t*  dz 


A  = 


voPg„Lv*P*cPs 


Where  vo  is  the  speed  of  sound  and  it  was  defined  as: 


(47) 

(48) 


(49) 


(50) 


(51) 


The  material  balance  in  dimensionless  form  was  defined  as: 

5r 

The  dimensionless  form  of  continuity  equation  was  written  as: 

dz  dZ; 


where: 


L 

The  dimensionless  pressure  is  defined  as: 

*  P 
P  = — 

Po 

Thus,  the  dimensionless  form  of  Ergun  equation,  Eq.40  was  defined  as: 


(52) 
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Eq.57  was  simplified  into  a  quadratic  form  written  as: 
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The  ideal  gas  law  in  dimensionless  fonn  was  defined  as: 


*  y 

p  =G 
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0  +  1 


where: 


(58) 
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(61) 
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Boundary  and  initial  conditions 

The  initial  conditions  are: 

r  =  0  0  <  £  <  1 :  0=  T°  ~T(:  ,  p* 

Tc 

The  boundary  conditions  are: 

r  >  0  £  =  0:  0  =  0,  p*=^- 

Po 

z  =  u  ^  =  o,  ^  =  o 


=  1,  rj  =  0,  v*  =  0 


Discrete  Crank-Nicolson  Implicit  Scheme  Method 

In  this  section  the  discrete  Crank-Nicolson  implicit  solutions  for  Eq.49  through  Eq.61 
were  formulated.  This  results  in  a  nonlinear  system  of  algebraic  equations.  A  clear  understanding 
of  these  nonlinear  algebraic  systems  of  equations  is  essential  for  describing  the  unique  way  in 
which  we  solve  them. 


Using  a  simple  first  order  in  time  (backward  Euler),  second  order  in  space  discretization  on  a 
uniform  grid,  a  discrete  version  of  energy  equation  (Eq.49)  was  written  as: 
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The  material  balance  (Eq.53)  after  discretization  gives: 
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Discretization  of  continuity  equation  (Eq.54)  gives: 
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The  discrete  Ergun  equation  (Eq.58)  was  written  as: 
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Discretization  of  ideal  gas  gives: 
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The  porosity  equation  (Eq.5.1 1)  was  written  as: 

C1  =  4>o  +  (70) 

Matrix-free  Approximation 

In  this  study  the  generalized  Crank-Nicolson  implicit  scheme  solutions  were  fonnulated 
into  matrices.  This  formulation  is  useful  in  solving  Equations  63  through  70.  After  simplifying 
Eq.63  gives: 


x;©”^1  +  s  •  ©r1  +  k;©;;/  -  r; 


38 


where: 
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Applying  the  boundary  conditions  to  Eq.71  result  in  a  matrix  (see  Fig. 38). 
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Figure  38.  Resulted  matrix  for  energy  equation. 


This  is  a  tridiagonal  system  i.e.  a  system  of  equations  with  finite  coefficient  only  on  the  main 
diagonal  (theZ/s),  the  lower  diagonal  (theX/s),  and  the  upper  diagonal  (theK/x). 


It  was  decided  that  the  velocity  of  the  gas  at  each  grid  point  and  time  r  +  A  r  has  to  be  evaluated 
using  Eq.68. 


*H+1 

V  ; 


4  a"  3 


(75) 


where: 
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As  a  result  of  this  mathematical  solution  to  velocity  of  the  gas,  the  continuity  equation  can  now 
be  used  to  solve  the  density  of  the  gas.  After  simplifying  Eq.66  it  becomes: 
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where: 
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Applying  the  boundary  conditions  for  i  =  1,2,... ,m  to  Eq.77  result  in  the  matrix  shown  in  Figure 
39. 
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Figure  39.  Continuity  equation  matrix. 

The  pressure  drop  along  the  cylindrical  sample  was  solved  using  the  ideal  gas  law  given  by 
Eq.69.  A  code  was  written  using  FORTRAN  program  to  perfonn  the  above  calculations.  The 
developed  code  was  incorporated  into  window  version  and  was  named  1-D  implicit  (closed 
system). 

Numerical  results 

Initially,  the  code  was  tested  for  non-reactive  system.  In  this  testing,  the  chamber  was 
pressurized  at  70  bars.  A  length  of  0.0254  m  of  the  cylindrical  sample  was  used.  The  porosity, 
initial  pressure,  particle  diameter,  viscosity  of  the  gas,  density  of  the  particle,  activation  energy, 
heat  of  reaction  per  mass  of  limiting  reactant,  thennal  conductivity,  and  mass  fraction  of  the 
limiting  reactant  were  set  to  0.4,  1.0  bar,  1.0e-5m,  2.29e-5  Pa.s,  2551.5  kg/m3,  227000  J/mol, 
5.85e6  J/kmol,  0.59  W/m.K,  and  0.60  respectively.  The  ignition  and  initial  temperature  were  set 
to  2000K  and  300K  respectively.  The  change  in  porosity  from  Eq.70  was  assumed  to  be  0.5.  In 
addition  to  the  mentioned  parameters,  the  heat  capacity  of  the  gas  during  the  testing  was  assumed 
to  be  constant  and  set  to  520.43  J/kg.K.  The  code  generates  pressure,  density,  and  temperature 
profiles  with  x  =  3.2  for  100  iterations  (see  Figures  40,  41,  and  42). 
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Figure  42.  Dimensionless  temperature  profile  for  non-reacting  system  with  At  =  0.005. 

Secondly,  a  reacting  sample  was  considered  in  which  the  pre-exponential  factor  was  set  to 
1.44el0s  .  The  pressure  in  the  chamber  was  reduced  to  50  bars  and  a  sample  length  of  0.0254  m 
was  considered.  The  sample  was  ignited  at  temperature  of  2000  K.  The  calculated  pressure, 
velocity,  density,  temperature,  and  the  fraction  of  reactants  reacted  are  shown  in  Figures  43-47. 


Figure  43.  Dimensionless  pressure  profile  with  At  =  0.005. 
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Dimensionles  velocity 


Figure  46.  Dimensionless  velocity  profile  for  50  bar  (725  psi)  with  At  =  0.005. 


Dimensionless  Length 


Figure  47.  Conversion  profile  for  50  bar  (725  psi)  with  At  =  0.005. 
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In  order  to  get  a  clear  understanding  of  this  system,  a  change  in  pressure  was  considered.  The 
sample  length  was  kept  at  0.0254  m.  The  pressure  inside  the  chamber  was  changed  from  70  bar 
(1000  psi)  to  140  bar  (2000  psi).  The  rest  of  parameters  were  kept  at  the  same  conditions.  The 
code  generated  the  pressure,  temperature,  density,  velocity,  and  conversion  dynamic  profdes. 
The  results  are  shown  on  Figure  48-57,  respectively.  The  velocity  ratio  was  plotted  at  various 
pressures  (725  -3000  psi)  and  the  results  are  shown  Figure  58. 


Figure  48.  Dimensionless  pressure  profile  with  At  =  0.004. 


Figure  49.  Dimensionless  temperature  profile  for  70  bars  (1000  psi)  with  At  =  0.004. 
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Figure  58.  Dimensionless  velocity  profile  at  various  pressures. 

As  can  be  seen  from  Fig. 35,  as  the  pressure  is  increased  from  725  psi  to  3000  psi  the 

dimensionless  velocity  increases  from  0.008  to  0.03  and  the  propagation  become  much  faster. 

However,  further  investigation  on  other  parameters  is  necessary  to  understand  this  system. 
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